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Abstract 

Quasi-stationary flows of gas accreting onto a compact center are analyzed in the framework of general- 
relativistic radiation hydrodynamics, under assumptions of spherical symmetry and thin gas approximation. 
Numerical investigation shows that luminosity, redshift and gas abundance are correlated. The gas can 
constitute up to one third of the total mass of brightest low-redshift sources, but its abundance goes down to 
1 /30 for sources with luminosities close to the Eddington limit. 
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We investigate a steady gas accretion onto a compact core in the framework of general relativity. 
The main goal of this letter is to show that bright sources — with the luminosity approaching 
the Eddington limit — must contain a significant fraction of gas. Our model assumes spherical 
symmetry, polytropic equation of state and thin gas approximation in the transport equation yj]. 

We use comoving coordinates t, r, < 6 < n, < (p < 2tt: time, coordinate radius and two 
angle variables, respectively. The metric is 

ds 2 = -N 2 dt 2 + Mr 2 + R 2 dQ 2 (1) 

where R denotes the area radius. The radial velocity of gas is given by U = 

The energy-momentum tensor reads = T^ v + T E where the baryonic part is given by 
T% - (P + pWfjUy + pg^ with the time-like and normalized four-velocity U M , U^U 11 = -1. 
The radiation part has only four nonzero components: T^ E = -p E = -T r r E and T E = T E r . A 
comoving observer would measure local mass densities, material p = T^U^Uy and radiation 
p E , respectively. The baryonic current is defined as f = poU^, where p is the baryonic mass 
density. Define as the unit normal to a centered (coordinate) sphere lying in the hypersurface 
t = const and k as the related mean curvature scalar, k = |V,n' = ~^d r R. The comoving radiation 
flux density reads j = U ll n y NTy E / = NT® E / V^. The baryonic matter satisfies the polytropic 
equation of state p = Kp T Q (with constants K and Y). The internal energy h and the rest and 
baryonic mass densities are related by p = p + h, where h = p/(T - 1). 
The equation 

V = (2) 

expresses the conservation of baryonic matter. 

There are four conservation equations resulting (due to the contracted Bianchi identities) from 
the Einstein equations, namely V^Ty 5 = -V M Ty E = F v (here v = 0, r). The quantity F v is the 
radiation force density and it describes interaction between baryons and radiation. The present 
formulation of general-relativistic radiation hydrodynamics agrees with that of Park y|], Miller 
and Rezzola [4] and (on a Schwarzschildean background) Thorne et. al [|50. 

One can solve formally the Einstein constraint equations G p o = S^T^o, arriving at ([2],[f]]) 



t= J,_^ + ( ;, (3 ) 



Above m(R) is the quasilocal mass given by 



m(R) = m-4nJ drr 2 (p + p E + -^j. (4) 



The integration in © extends from R to the outer boundary 7?oo. A ball of gas is comprised between 
a hard core of a radius Rq and sphere Soo of a radius R^. Its external boundary is connected to 
the Schwarzschild vacuum spacetime by a transient zone of a negligible (due to special transitory 
data) mass. Thence the asymptotic mass M is approximately equal to m{R oa ). Similar picture 
emerges in the recent construction of quasistars |7|]. 

In an alternative, polar gauge foliation, one has a new time ts(t, r) with d ts = d t - NU8r. The 
expression 4-nNkR 2 + ( j ) j + 2Up E /k^ represents the radiation flux measured by an observer 
located at R in coordinates (t s ,R). One can show that 

d ts m{R) = -4n(NkR 2 + ( j) ) + ) + NUR 2 (p + p)j . (5) 

The mass contained in the annulus (R, Rco) changes if the fluxes on the right hand side, one directed 
outward and the other inward, do not cancel. 

The local baryonic flux will be denoted as M = -AnUR 2 p G and its boundary value reads M M . 
The accretion process is said to be quasistationary if all relevant observables measured at R are 
approximately constant during time intervals much smaller than the runaway instability time scale 
T = M/Moo. Analytically, we assume that d, s X = (d t - NUd R )X = for X = p , p, j, U. . . 

The above assumptions imply, in the thin gas approximation [1], that F = and the radiation 
force density has only one nonzero component F r = KkNp j. The only direct interaction between 
baryons and radiation is through elastic Thomson scattering, k is a material constant, depending 
in particular on the Thomson cross section cr, k = cr/ (4-7mi p cy c is the speed of light and m p is the 
proton mass. 



The full system of equations in a form suitable for numerics has been derived elsewhere 111 ill . 
It consists of: 

i) The total energy conservation 

r-i . P E , / u 2 \ 

MN + 2MN— = AnR 2 jNk 1 + — + C; (6) 

r - 1 - a 1 po \ k J 

The constant C is the asymptotic energy flux flowing through the sphere of a radius i?oo (see 
©). 

ii) The local radiation energy conservation (below a = is the speed of sound) 

2m(R)\ N d i , FX / F ,dU 



R 



2k (jU - kp E ) ^ + SnNRk (/ - jp £ |) . (7) 
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iii) The relativistic Euler equation 



d , 2 
— In a 

dR 



r - 1 - a 2 

a 2 _ Ml 

a k 2 



x 



r- 1 



iv) The baryonic mass conservation 

dU 

1r 



U d , , 2U AnR 

— lna 2 + -. 

r - 1 - a 2 dR R k 



(8) 



(9) 



v) The equation for the lapse 



dN m / T- 1 -a 2 d , / , ,x\ 
— = N\Kj + — lnfr- 1 - a 2 ) . 



(10) 



Equations © — [TOl) constitute, with and m(7?) given by © and ©, the complete model used in 
numerical calculations. The asymptotic data are such that a 2 ^ » M/R^ » f/^,, which guarantees 
the fulfillment of the Jeans criterion for the stability (see a discussion in [80 and studies of stability 
of accreting flows in newtonian hydrodynamics [90), suggesting in turn the stability of solutions. 
One can put _/'«, = pf,. The total luminosity is well approximated by L = AnR^joo and it should 
be related to the accretion rate by the formula 
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N(R Q ) 
k(R ) 



2m (R ) 
Ro 



M 



(ID 



The last formula is justified by two arguments, i) In the nonrelativistic limit one gets a = \<p (R ) |, 
where is the newtonian potential. Eq. (fTTT) states now that the binding energy is transformed 
into radiation with the implicit assumption that the heat capacity of the core is negligible, ii) The 
condition of stationarity implies the existence of the approximate time-like Killing vector and it 
appears that a gives the standard measure of the gravitational redshift. If stationary observers 
detect a>o at Rq and to at infinity, and 1 /to «: 2M (the geometric optics condition — see for 
a discussion) then to = (1 - a) to . Thus a can be regarded as a proper binding energy and again 
one arrives at formula (fTTT) . 
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Let us remark that from above definitions and the equation ([8]) one infers Lo < 4nM/K, for 
accretion solutions; it is notable here that the limiting luminosity involves the total mass M instead 
of the mass m(R ) of the central core. 

The triple of independent boundary data can consist of a, L and a 1 ^. These quantities can be 
determined from observations of highest redshift, total luminosity and asymptotic temperature, 
respectively. Then one chooses _/«, = pf, = L$l (4nR 2 J}, and the mass accretion rate M = L /a. 



These data specify supersonic flows up to, possibly, a bifurcation 111 111 . In the case of subsonic 
flows another boundary condition is needed, for instance the asymptotic baryonic mass density 

Poo- 

Eqs. (TTlfTOl) are in the evolution form. Numerical calculation starts from the outer boundary R^, 
taking into account Eq. ©, and evolves inwards until the equality a = 1 - Jl - is met 
at some R, denoted as R and being regarded as the radius of the compact core. The numerical 
integration employs the 8th order Runge-Kutta method Jl2]. Choosing p 0oo at random one either 
obtains no solution at all or a subsonic solution. Using the bisection method and automating the 
search process, one can obtain a boundary of the solution set (on the plane L - pooo), which 
(interestingly enough) appears to bifurcate from a brightest flow. This boundary will be called 
later on as the bifurcation curve. 

We choose M /M = 5.95496 x 10" 7 , where M is the Solar mass. In the standard gravitational 
units G = c = 1 and in the scaling M = I one gets k = 2.1326762 x 10 21 (M /M), that is 
k = 1 .27 x 10 15 . The size of the system is Roo = 10 6 . The speed of sound is given in successive runs 
bya^ = 4xl(T 3 ,4xl(r 4 ,4xl(r 5 . The Eddington luminosity reads L E = 4nM/x = 9.9847x 10~ 15 . 

Figures |T3-[3] show accreting solutions on the luminosity-(mass of the central core) diagram 
for a = 25 x 10~ 4 , 0.5, 0.9, respectively. Each figure depicts solution sets for three different 
values of the asymptotic speed of sound a^. For small a there can exist two accreting solutions 
possessing sonic points, with asymptotic densities p 0oo i and pooo2- Subsonic flows then exist for 
each poco e (poooi,Pooo2)- Subsonic solutions are not specified uniquely for given boundary data 
but the length of the interval of allowed values of the asymptotic baryonic density p 0oo becomes 
shorter with the increase of L . 

For larger L and/or a the bifurcation curve can consist either of subsonic or supersonic flows 
and its interior consists of subsonic solutions [11]. The brightest system coincides, as before, with 
the bifurcation point and it is unique. The luminosity of the bifurcation point increases with the 
decrease of the asymptotic speed of sound and it goes up with the increase of a. Its gas abundance 
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FIG. 1: Small binding energy, a = 0.0025. Three asymptotic values of the speed of sound, ai = 
0.004, 0.0004, 0.00004. Two bifurcation branches encompass the set of subsonic flows. The abscissa shows 
the luminosity and the ordinate shows the mass of the compact core. 

depends both on luminosity and redshift. 

The gas abundance for a = 25 x 10~ 4 decreases from almost 1/3 at ai = 0.004 to 1/5 at 
ai = 0.00004, as shown on Fig. 1. This value of a implies 2M(R )/R « 0.005. Interestingly, the 
abundance 1/3 can be shown analytically to characterize those general relativistic polytropic flows 



without radiation that maximize the accretion rate LI 3] and low luminosity newtonian sources 



[|8|]. The case of a = 0.5 corresponds to a very compact central object with 2M(R )/R « 0.75, 



close to the Buchdahl limit H14H . Fig. 2 demonstrates that gas contribution equals about 0.16 for 
ai = 0.004 and goes down to 0.04 for ai = 0.00004. When a = 0.9, then 2M(R ) /R Q ~ 0.99 at the 
surface of the compact central object, well beyond the Buchdahl limit. Only exotic matter violating 



standard energy conditions can be responsible for such compact bodies II 1511 . We find from Fig. 3 



that gas contribution to the mass changes from 1/8 (ai = 0.004) to 1/30 (ai = 0.00004). 
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FIG. 2: a = 0.5, a* - 0.004, 0.0004, 0.00004. The axes are as in Fig.l. 

It is clear that a reformulation of the problem would allow one to estimate the mass of an 
isolated system, assuming that the mass of the central core is known. This can be possibly applied 
to Thorne-Zytkow stars. 

One can show that in models with test fluids the gas density is bounded from below. In partic- 
ular, in a Shakura model and for general relativistic systems with low luminosity and redshift, the 



bound is provided by a supersonic flow [11]. The full general-relativistic analysis reveals a new 
qualitative effect. Namely in steady (sub-or supersonic) accretion solutions the gas abundance is 
bounded both from below and from above by bounds that depend on the redshift and luminosity. 
This is a clear demonstration of the importance of backreaction in accretion processes. 

The triple of observables a, L , does not specify accretion completely, but for high lumi- 
nosities the remaining freedom (in choosing the asymptotic baryonic mass density) is severely 
restricted and the brightest flow is unique. It is interesting to note that the concept of Eddington 
luminosity still applies - in the light of our data - in the general relativistic case. In all stud- 
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FIG. 3: High binding energy, a = 0.9, - 0.004, 0.0004, 0.00004. The axes are as in Fig.l. 

ied examples we have L < L E ; this is supported also by an analytic argument, discussed above. 
Now L E = AtiMIk; it is the global mass rather than the mass of the compact core, that enters the 
expression for the Eddington luminosity. 

Numerical data show that gas can be abundant in quasi- stationary accreting systems. Brightest 
systems can possess even 33% of gas for small redshifts and still more than 10% of gas for a = 0.9. 
It is an open problem whether this conclusion is true in nonspherical steady flows. 
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